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ABSTRACT 

In this paper we show that the high energy 7-ray flux in the GeV domain from 
mature pulsar wind nebulae (PWN) scales as the change in rotational kinetic en- 
ergy — Q 2 )/2 since birth, rather than the present day spindown power IQQ. 
This finding holds as long as the lifetime of inverse Compton emitting electrons 
exceeds the age of the system. For a typical 7 -2 electron spectrum, the predicted 
flux depends mostly on the pulsar birth period, conversion efficiency of spindown 
power to relativistic electrons and distance to the PWN, so that first order es- 
timates of the birth period can be assessed from GLAST/LAT observations of 
PWN. For this purpose we derive an analytical expression. The associated ("un- 
cooled" ) photon spectral index in the GeV domain is expected to cluster around 
~ 1.5, which is bounded at low energies by an intrinsic spectral break, and at 
higher energies by a second spectral break where the photon index steepens to 
~ 2 due to radiation losses. However, if the spectral parameters deviate from the 
abovementioned assumptions, we can combine GLAST/LAT observations with 
multiwavelength data from radio to TeV 7-ray observations to improve such birth 
period estimates. From the ATNF pulsar catalog we expect a total of 50 galactic 
PWN to be detectable if the ratio of birth period to current period is < 0.5, but 
30 for a ratio < 0.75. Mature PWN are expected to have expanded to sizes larger 
than currently known PWN, resulting in relatively low magnetic energy densities 
and hence survival of GeV inverse Compton emitting electrons. Whereas such 
a PWN may be radio and X-ray quiet in synchrotron radiation, it may still be 
detectable as a GLAST/LAT source as a result of the relic electrons in the PWN. 

Subject headings: pulsars: general — Radiation mechanisms: non-thermal - 
ISM: jets and outflows — ISM: supernova remnants — gamma rays: theory 
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Introduction 



The most reliable method of determining the initial spin period of a pulsar is to 
have a historical association of the s upernova explosion, combined with a m easurement 
of the pulsar braking index: For Crab, iLyne. Pritchard & Graham-Smithl (119931 ) estimated 
a birth period of Pq = 19 ms, whereas a value of Pq = 63 ms was determined for PSR 



B 15 09-58. The pulsar curr e nt tec hnique was introduced by iPhinnev fc Bland ford (1981) 



and Viyekanand & Narayanl (1981), but suffers from small number statistics ( Lorimer et al. 



19931 ). iNarayanl (119871 ) and lEmmering fc Chevalierl (119891 ) proposed a technique based on 
the pulsar luminosity function technique, injection and the presence of pulsar wind nebulae 
(PWN), so that injection at Pq ~ few hundred milliseconds would result in the absence of 



a PWN. Ivan der Swaluw fc Wul (120011 ) have shown that the pulsar wind nebula radius is 
a function of the birth period of a pulsar, and that the ratio of this radius relative to the 
supernova remnant forward shock radius provides a measure of the birth period as well. 

In this paper we show that high energy gamma-ray observations of pulsar wind nebulae 
in the GeV domain can also give us a measurent of the birth period of the associated pulsar 
wind nebula. The reasoning behind this is simple: the relatively low energy leptons radiating 
GeV 7-rays via inverse Compton scattering on the cosmic microwave background (CMBR) 
and the galactic radiation fields do not suffer from radiation losses as e.g. X-ray synchrotron 
emitting electrons do. In this case the GeV flux should be proportional to the total energy 
in leptons ejected since birth, which in turn is a fraction rj of the total rotational kinetic 
energy E rot = ~~ f2 2 )/2 released since birth. Here Qq — ^tt/Pq and Q are respectively 
the birth and current pulsar angular frequencies. 

Normally such a technique would suffer from uncertainties in PWN parameters, but we 
will show that if we assume a particle injection spectral index of 2 above the radio domain, 
the expected GeV flux can be expressed as a simple analytic formula, in which case the 
uncertainty in P scales linearly with the uncertainty in the distance, but less sensitively 
on other less known parameters such as rj, the maximum particle energy 7 max ' m eC 2 and 
the intrinsic spectral break energy 7f,m e c 2 above the radio domain. The latter typically 
divides the injection spectrum into a ~ 7~ 2 component above this break, but a much harder 
spectrum at 7 <C 7& to account for the flat spectra seen from radio PWN. The requirement 
for such an intrinsic break in the Crab nebula and other PWN was discussed respectively by 
Atoyan fe Aharonianl (119961 ) and I Chevalierl (120041 ) . 



Whereas such an analytical expression helps us with population studies and to get ap- 
proximate estimates for Pq, a GLAST/LAT detection with a clear multiwavelength PWN 
association should not rely on this expression to obtain an accurate estimate of Pq: The 
multiwavelength synchro- Compton spectrum should give us more information about parti- 
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cle spectra, break energies and conversion efficiencies through proper numerical modeling, 
although the principle remains the same as expressed below. 



2. The electron spectrum contributing to the GLAST/LAT spectral range 



We adopt the injection spectrum of IVenter fc de Jagerl (120061 ) for electrons at the pulsar 
wind shock 

t\ = ( Qo(t)(lhb)- pi for 7 < lb \ 
W ' tj ^ Qo(t)(l/l b r P2 for lb < 7 < 7max J ' 1 } 

with 7;,m e c 2 the intrinsic break energy. For p\ ~ 1 we retrieve the typical flat radio spectra 
of PWN where F v oc v~ a with a ~ 0. The high energy spectrum is steeper and can extend 
up to ultra high energies (7max^ e c 2 ~ 10 14 eV to > 10 15 eV). The corresponding spectral 
index p 2 ~ 2 as a result of particle acceleration at or near the pulsar wind shock. This is 
supported by observations at X-ray and very high energy 7-ray energies: At the pulsar wind 
shock (radius R s ) we observe X-ray and VHE 7-ray photon indices T s = (p 2 + l)/2 ~ 1.5, 
whereas for r 3> R s the electron spectrum cools due to synchrotron l osses (p 2 i ncrea ses by 



unity), resulting in a value of T x = (p 2 + 2)/2 ~ 2 as illustrated by Ide Jagerl (120081 ). who 
combined the results from several PWN in the form of a plot of T x vs T s . This reveals a 
cluster near T s ~ 1.5, evolving to T x ~ 2. In this paper we will assume p 2 = 2, resulting in 
a logarithmic suppression of uncertainties in % and 7 max , which is important for a relative 
accurate determination of the birth period. Note however that there are deviations: There 
are a few cases for which p 2 ~ 3, in which case a multiwavelength/numerical modelling 
approach is called for, as discussed in the Introduction. 

In general we find a critical energy 7 r ad^eC 2 , such that radiation losses are unimportant 
for 7 <C 7rad (the "uncooled" domain) and R s < r < -Rpwn (with -Rpwn the PWN radius), 
in which case the total observed spectrum is reflected by the injection indices p\ and p 2 - For 
7 3> 7rad (the "cooled" domain) , the corresponding electron spectral index steepens to p 2 + 1 
due to radiative cooling, in which case total nebular photon spectra of the form E~ 2 are 
observed. An estimate of 7 ra d can be obtained by setting the timescale for radiative losses 
(assuming isotropic lepton pitch angles) equal to the age T age of the system and rewriting 
this electron energy in terms of its inverse Compton scattered 7-ray energy giving (assuming 
scattering on the CMBR): 



£; 7 (rad) = -(2.7kT)7 r 2 ad -(0.50TeV)( 777 - 7 ,) ( ) . (2) 



-(2.7kT) 7 2 ~ (0.50 TeV) ( — ■ , , 
3 V Jhad K 'yiO^Gj ViOkyr, 

For B < 3/xG this expression must be corrected for inverse Compton losses on the CMBR. 
The spectral breaks of the two mature H.E.S.S. sources HESS J1825-137 (i? 7 (rad) ~ 2 TeV; 
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'agc 



20 kyr - lAharonian et al 



2006bh ) and the cocoon of Vela X (£ 7 (rad) ~ 13 TeV; 



r agc = 11 kyr - lAharonian et al.l (j2006al )) can be explained with field strengths between 4 



and 5/iG. Furthermore, the relati vely faint X-ray synchrotron flux from these sources support 
such low field strengths (see also Ide Jagerl tooj)). The newly det ected PWN HESS J1718- 
385 also shows a spectral break at 7 TeV (lAharonian et al.l 120071 ) and the discovery of an 
X-ray PWN around the associated pulsar PSR J1718-3825 implies a field strength < 5/iG 
( iHinton et al.ll2007l ). For a characteristic age of 90 kyr the predicted break energy is however 
in the multiple GeV range. If the birth peri od is not much less t han the current 75 ms period, 
it would be possible to reproduce the break (IHinton et al.ll2007l ). but for this we need to know 
the pulsar braking index. 

Assuming evolved PWN with radiation breaks not lower than 10 GeV and B ~ 5/iG, 
GLAST/LAT should be able to detect PWN between the intrinsic break energy and ~ 10 
GeV with a ("uncooled") photon index near 1.5 if the age 



r age < (280 kyr) 



B 



5/iG 



(3) 



Having set the constraints for detecting the specified injecti on spectral indices, we can now 
proceed with the high e nergy flux calculations: Following iSefako fc de Jagerl (120031 ) and 
Venter &: de Jagerl (120061 ). the energy equation can be written in terms of the present day 
spindown power L(t) giving 



m e c 2 J Q(7, t)jdj = rjL(t). 
From this expression the normalisation constant Qo(t) for p 2 = 2 follows: 



(4) 



Qo{t) 



yL(t) 



with a 



2 -Pi 



+ lu- 



ll 



(5) 



and rj the conversion efficiency of spindown power into energetic particles. Since we observe 
at electron energies well below 7 ra( j (assuming escape losses are also unimportant), the total 
nebular particle spectrum at r > R s can be derived from the transport equation to give 



rfiV(7,T, 



agcy 



G?7 



QVY,t)dt — ^ 2 7 . 



am e c^ 



(6) 



The time integral J L(t)dt up to the present age is given by the change in rotational kinetic 
energy AKE rot = /(fig - fi 2 )/2 = 2tt 2 /(1/P 2 - 1/P 2 ). 
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3. The expected 7-ray flux and the sensitivity of LAT to birth periods 

The differential gamma-ray spectru m resulting from the inverse Compton scattering of 



the CMBR (with T = 2.7K) is given bv Blumenthal fc Gouldl (11970 ). We then replace the 



electron spectral normalisation constant with the abovementioned rotational kinetic energy 
to give a total differential photon rate of 



dN~ ?~o i] 



dtdE^ 7r/i 3 m e c 4 a 



^AKE rot (kT) 7 ^F(p 2 )E- 1 - 5 . (7) 



Note the photon spectral index of 1.5 resulting from el ectrons radiati n g wit h a p 2 



spectral index. The function F(p 2 ) = 5.3 for p 2 = 2. Ide Jager et al.l (119951 ) introduced 
the necessity of including the contribution of far infrared (FIR) photons from galactic dust 
grains (at a temperature of ~ 25K) to the inverse Compton process for PWN for which 
external Compton scattering is more important than the synchrotron-self-Compton process. 
In general, if Ui (in units of eV/cm 3 ) is the energy density of an ambient target photon field 
with corresponding temperature T,, the additional flux in the Thomson limit due to species 
i wo uld then be given by ( i n the energy range where p 2 = 2 contribute to the 7-ray flux - 
from iBlumenthal fc Gouldl Jl97oh ) 



dNjW) ^ g 5T - 1/2u dN,(T = 2.7K) 
dtdE 1 ' ' dtdEry 



With this derivation we assumed that the energy density at temperature T, is a constant 
times that of a pure black body at the same temperature. The energy density of this FIR 
is U ~ 1 eV/cm~ 3 within a galac tocentric radius of ~ 4 kpc, decreasing to ~ 0.3 eV/cm -3 
towards the solar neighbourhood ( Porter et al.ll2006l ). From this expression for the 7-ray flux 



it is clear that high temperature photon fields do not contribute, unless the corresponding 
photon energy density is relatively large. We will therefore only add the contributions for 
the CMBR and FIR dust fields. 

By selecting p 2 = 2, the constant a in Qo(t) can vary between 6 and 12 for pi ~ 1, 
7 max mc 2 = 10 14 eV (Vela-like) to 2 x 10 15 eV (Crab-like) and a break energy of 7f,m e c 2 
between 10 11 eV and 3 x 10 12 eV to reproduce the observed break frequencies near 10 13 
Hz and 10 16 Hz respectively in the Crab nebula. We will normalise a to a value of 10 by 
setting a = 10ai . In general, since the birth period scales as a _1//2 , uncertainties in the 
abovementioned parameters are significantly suppressed since we take the square root of a 
logarithm. 

We can now write the total high energy 7-ray spectrum from a PWN in the GLAST/LAT 
GeV domain (i.e. where 7b < 7 < 7 ra d), where the electron spectral index over the whole 
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nebula remains ~ 2, so that the photon index is ~ 1.5. If the pulsar is also sufficiently spun 
down (P Q < P), we can neglect the P-term in AKE rot , so that 

dN 1 _,„,„-s(l + l-3f») r =C ,_ L5 



^ 7 - 4xl0 " 8 ^cfw 5 < 9 > 

in units of cm~ 2 s _1 GeV . Here I45 is the moment of inertia normalised to 10 45 g.cm 2 , 
w hereas the birth period is nor malised to 40 ms: Pq = (40ms)P4o, following the finding 
of Ivan der Swaluw Wul (120011 ) that many pulsars in composite remnants were born with 
Pq ~ 40 ms. From the integral flux above 1 GeV we can calculate the LAT sensitivity to 
detect a given birth period: For a 50 hr, 5a flux sensitivity of 2 x 10 -10 cm _2 s~ 1 above 1 
GeV it is clear that LAT should be sensitive to birth periods in the range 

P <(44 ms )( 1 + , 3fe ,^(^) 1/2 (_i_)- 1 . (10 ) 

For example, for a source at a distance of 10 kpc, but within a galactocentric radius of 4 
kpc where U25 ~ 1 eV/cm~ 3 , GLAST should be able to detect birth periods less than 67 
ms, whereas for a source at a distance of 2 kpc (U25 ~ 0.3 eV/cm~ 3 ) birth periods as long 
as 260 ms can be detected, assuming a realistic conversion efficiency of 77 = 0.3. 

Assuming birth periods Po < eP (with e < 1), we can use this flux prediction (in- 
cluding the P-term) to determine the number of pulsars with expected GeV fluxes above 



the G LAST /LAT sensitivty threshold. Using the ATNF pulsar catalog of Manchester et al. 



(120051 ) . we select all pulsars with characteristic ages less than 280 kyr (according to the con- 
dition set by Eqn. (j3J)) and for each of these pulsars we calculate the integral flux above 1 
GeV from Eqn ([§]) for a given period and selected value of e (the current period was included 
for this calculation). Note that we only take the total number from this catalog and did 
not correct for selection effects such as dispersion measure constraints and off-beam pulsars. 
Adding such potentially unseen pulsars should increase the numbers presented below. The 
result from the existing catalog is a total of 50 detectable pulsars for e = 0.5, but only 
30 for e = 0.75. Note however that the true birth periods are not expected to scale as a 
fixed fraction e of the current period - the only purpose of these estimates is to show that 
a relatively large number of PWN are expected to be visible even if P differ by only 25% 
from P. 



4. Discussion and conclusion 



In this paper we derived the expected GeV 7-ray spectrum of PWN under the assump- 
tion that the particle spectrum in the uncooled domain (7 < 7 md ) is p 2 ~ 2, as observed 
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from the injection spectra at most PWN shocks. When observing mature PWN (i.e. those 
for which the age is significantly larger than the spindown timescale), we arrive at a fresh 
new perspective w.r.t. to the different multiwavelength domains: 

For the GLAST/LAT domain above 1 GeV, the observed 7-ray luminosity would scale 
as the change in rotational kinetic energy I(Qq — Q 2 )/2, rather than the present day spindown 
power IQQ. This is because the electrons contributing to the LAT energy domain survive 
since the earliest epoch when the spin period was close to the birth period. We therefore 
predict that the GeV luminosity and IQQ would be uncorrelated. There is however a caviat 
to this argument: The field strength during earlier epochs is expected to be larger, so that 
GeV emitting electrons may have been burnt off during an earlier epoch of h igh B. However, 



for the PWN of PSRB15 09-58 the n ebular field strength is between 8/xG (Idu Plessis et al. 



19951 ; iGaensler et al.ll2002l ) and 17/^G ( lAharonian et al.ll2005l ). given an age of ~ 1, 700 years, 
whereas Eqn (jSJ) requires a field strength of B < 64/xG after 1,700 years to force a radiation 
spectral break above 10 GeV. It is thus clear that this requirement is easily achieved for this 
example. The same can be shown for other relatively young VHE emitting PWN, such as 
Kes 75, G21. 5-0.9 and GO.9+0.1, except for the Crab Nebula for which B~2x 10~ 4 G after 
950 years. In general, the expansion history of the PWN, given the environment into which 
the associated supernova remnant forward shock expands, is expected t o determine the time 
evolution of PWN field strength. See also iReynolds fc Chevalier! (I1984I ). 



For the very high energy 7-ray (ground based) domain, we detect i nverse Compton 
radiation from electrons around the radiation maximum 7 ra( j as reviewed by Ide Jagerl (120081 ) 
for PWN older than 10 kyr (see the disc ussion on Vela X and HESS J1825-137). This is 
possibly also the case for HESS J1718-385 (lAharonian et al.ll2007l : lHinton et al.ll2007l ). Below 
this maximum the photon spectral index is near 1.5 due to relatively old electrons ejected 
since birth. This hard component should extend into the GLAST/LAT range as discussed 
in this paper. However, electrons with energies 7 3> 7 ra d have lifetimes shorter than the 
age of the pulsar, so that they are expected to reflect the current spindown power. Thus, 
for the VHE domain some correlation with IQQ may be expected, with the strength of the 
correlation depending on the magnitude of 7 ra d- 

In the X-ray domain we observe only synchrotron radiation from freshly ejected electrons 
(i.e. 7 3> 7rad), in which case we exp e ct a relatively tight correlation with the current 
spindown power. In fact, ICheng et al.l (120041 ) isolated the PWN X-ray contribution and 
found a correlation between X-rays and spindown power of the form Lx oc (IQQ) 1A±0A . 



Recently iHuang fc Wul (120031 ) derived birth periods for normal pulsars in the range 
Pq ~ 0.6 to 2.6 ms. For such birth periods this study would predict very bright PWN at 
a distance of 10 kpc - even by EGRET standards, unless t) <C 1 during early epochs after 
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birth. However, observations of the Crab nebula does not support a low value of ry at birth. 

Mature (Vela-like) PWN appear to have field strengths below 5/xG. If we also require the 
spectral turnover for such PWN due to radiative losses to be above 10 GeV, the corresponding 
age of the PWN should be < 300 kyr as derived from Eqn. [3j Assuming that mature PWN 
evolve to such low field strengths in general, we infer that ~ 50 PWN should be detectable 
for GLAST/LAT if we assume a birth period Pq < 0.5P, whereas this number reduces to 
30 for P < 0.75P. More may be visible if a significant amount of pulsars are off-beam, 
in which case we may not see any pulsed counterpart. Another implication of an evolved 
PWN with age near 300 kyr is that the current spindown power and magnetic field strength 
are too low for a detectable synchrotron nebula, whereas the relatively low energy relic 
lepton component would still inverse Compton scatter external radiation fields into the GeV 
domain, which would lead to a population of "dark high energy 7-ray sources", i.e. those 
without multiwavelength counterparts. 

To avoid source confusion and claim sources with confidence it would be important 
to follow the LAT spectrum into the VHE 7-ray domain. The advantages for such an 
association would be the following: (a) We detect the spectral turnover associated with 7 rad , 
(b) the positional agreement of the images between GeV and and the higher resolution TeV 
energies makes the claim more significant, and (c) the size is expected to shrink for photon 
energies 3> 3kTj^ d , with the i mage centroid convergin g towards the associated pulsar as 



observed from HESS J1825-137 (lAharonian et al.l l2006bl ) . Finally, to improve the accuracy 



in determining P Q , one should also add multiwavelength spectral information, if observable. 



The author would like to thank Mallory Roberts, Arache Djannati-Atai and several 
members of the H.E.S.S. working group for supernova remnants, pulsars and plerions for 
useful discussions on the VHE 7-ray and multiwavelength properties of plerions. 
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